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1 Introduction 

The role of self-decomposable (SD) distributions in first order autoregressive 
(AR(1)) models of the form 

X n = cX n _ 1 + e n , (1.1) 

described by random variables (r.v.s) {X n ,n e Z}, innovations (i.i.d. r.v.s) 
{e n } and c G (0, 1) such that for each n, e n is independent of X n _i, has been 
discussed by many authors, see e.g. Bouzar and Satheesh (2008) and the refer- 
ences therein. Recently Kozubowski and Podgorski (2010) has introduced the 
notion of random self-decomposability of distributions on the reals motivated 
by stationary solutions to the AR(1) model 
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e n , with probability p, 

-X-n — \ U-2J 

cX n _i + e n , with probability (1 — p), 

described by r.v.s {X n ,n e 2"}, innovations {e Tt } and c e [0, 1] such that for 
each n, e n is independent of X n _x- 

Definition 1.1 A charcteristic function (CF), of a probability distribution, 
ip{t) is randomly self-decomposable (RSD) if for each p,c G [0, 1] there exists 
a CF ip c ,p(t) such that 

VK*)=^*(*){p+ (!-*#(<*)}. (1.3) 

Kozubowski and Podgorski (2010) then discusses the relation of RSD laws 
to SD laws and geometrically infinitely divisible (GID) laws. In proposition 
2.3 they prove, in an elegent manner, that the class of RSD laws equals the 
intersection of the classes of GID laws and SD laws. They also discuss a variety 
of examples. We need the following in our discussion. 

Definition 1.2 Harris(a, k) distribution on {1, 1 + k, 1 + 2k, ....} is described 
by its probability generating function (PGF) 

P(s) = t -, ~t k > integer and a > 1. (1.4) 

w {a - (a - l)s k y/ k v ' 

Definition 1.3 A CF ij}{t) is Harris-ID (HID) if for each p e (0, 1) there exists 
a CF ip p (t) such that 

Theorem 1.1 (Satheesh et al. (2008)) A CF ip(t) is HID iff 

*<" = a-iogW/' (L6) 

where > integer and /i(t) is some CF that is ID. 

When k — 1 Harris distribution becomes the geometric(])) distribution on 
{1,2,...} withp= i. For more on this distribution see Sandhya et al. (2008). 
Certain aspects of HID laws and generalized AR(1) models have been discussed 
in Satheesh et al. (2008). In section 2, the notion of RSD is generalized, 
its relation to SD laws and HID laws are presented and its connection to a 
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stationary generalized AR(1) model is given. The notion is then extended to 
Z + -valued distributions in section 3. We closely follow the development in 
Kozubowski and Podgorski (2010). 



2 Generalizing RSD distributions 

Remark 2.1 In the paragraph after their Proposition 3.1 Kozubowski and 
Podgorski (2010) state that AR(1) processes described by (1.2) cannot be 
constructed with either (general) gamma or Gaussian distributions for X n as 
neither of them are GID although both are SD. However, it should be noted 
that gamma (a, A) distributions (equation (2.10)) are GID if a < 1, see e.g. 
Yannaros (1988) or Sandhya (1991). 

Definition 2.1 A CF ip(t) is Harris-RSD (HRSD) if for each c G (0, 1] and 
each p G [0, 1) there exists a distribution with CF ip c , P if) such that 

m = ^ P (t){ P + (i - p)^ k {ct)y k . (2.i) 

Remark 2.2 With the above nomenclature the RSD defined by Kozubowski 
and Podgorski (2010) is geometric RSD (GRSD) because it bridges the notions 
of SD and GID where as our definition bridges the notions of SD and HID. 

When p = equation (2.1) reduces to 

i;(t) = ^{ct)^j c {t) (2.2) 

where tl> c (t) = ip c ,o(f), that is ip{t) is SD. On the other hand when c = 1 
equation (2.1) becomes 

^{t) = ^ p {t){p+{l-p)^\t)}^ k . (2.3) 

where ip p (t) = ipi, p (t). Solving for t/j(t) we get 

m = i 7 ^Wvrr7r ;a = -- ( 2 - 4 ) 

VK) {a-(a-l)^(t)}V*' p 1 ; 

That is if>(t) is HID. 

Denoting the classes of HRSD, SD and HID distributions by C HRS d, Csd 
and Chid the above discussion shows that Chrsd C Csd H Chid- In the next 
Proposition we show that we have equality here. 
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Proposition 2.1 We have Chrsd = Csd H Chid- Further, whenever the CF 
ijj(t) £ C H rsd, the CF ip CiP (t) in (2.1) can be written as 

^M*) = i) c {t)4 P {ct) (2.5) 
where ^c(^) and are given by 

<M*> = { p +(i-p)V'(«)}v= (2 ' 7) 

Proof. If the CF ^(t) is SD then for each c £ (0, 1] the function ip c (i) in (2.6) is 
a genuine CF and similarly if vp(t) is HID then for each p £ [0, 1) the function 
^/>p(£) in (2.7) also is a genuine CF. Consequently (2.5) is a well defined CF 
and hence (2.1) holds, proving the assertion. 

Now let us consider a generalization of the AR(1) sequence (1.2). Here 
{X n } is composed of k independent AR(1) sequences z = 1, 2, . . . k and 

where for each n, {Y ni } are independent. That is, for each n, X n = Ei=i Y n ,i 
and e n = Ei=i e n,i where is an i.i.d sequence and similarly {e n ,i} is also 

an i.i.d sequence, k being a fixed positive integer. Further, it is also assumed 
that for each n, e Uji is independent of Y n -\^ for all i — 1,2, .. . k. Situations 
where such a model can be useful have been discussed in Satheesh et al. (2008). 

_ ( Ef=i e «,i, with probability p, ^ ^ 

i=i [ELi cY n-i,i + E!=i e n,i, with probability(l - p). 

In terms of CFs and assuming stationarity we get 

Mt) = A(t){p+(l- P )^ k Y (ct)Y /k . (2.9) 

The following Proposition is now clear. 

Proposition 2.2a If {Y ni } describes the AR(1) model (2.8) that is stationary 
for each c £ (0, 1] and p £ [0, 1) then the distribution of {Y nyi } is HRSD. 

Now proceeding as in the proof of Theorem 2.2 in Bouzar and Satheesh 
(2008) we have the following converse to Proposition 2.2a. 

Proposition 2.26 If ^y{t) is a CF that is HRSD with ip c ,p(t) = i>e(t) for each 
c £ (0,1] and p £ [0, 1) then there exists a stationary AR(1) model described 
by (2.8) with the CF of {Y n i } and ip e (t) that of the innovations {e n ^}. 
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Example 2.1 Gamma(|,A) distributions has CF 

m = (1 _* At )iA ' k > inte s er ' A > ( 2 - 10 ) 

is HID. Further since it is also SD, this distribution is HRSD. 

Example 2.2 Let the CF ip(t) be Harris-sum-stable for every c G (0, 1). Then 

^'^■ (.-(.-iVw)-" ' (2 ' n) 

The second factor on the RHS is also a genuine CF being a Harris-sum of ip(t) 
where the support of this Harris distribution is Z+. Thus ip(t) is SD. Hence 
if we take h(t) as a stable CF in Theorem 1.1 then ip(t) in (1.6) is Harris- 
sum-stable for every c G (0,1) and we have a general procedure to construct 
CFs that are HRSD. With k — 1 above, we have the corresponding geometric- 
sum-stable laws and a procedure to construct the examples in Kozubowski and 
Podgorski (2010). 

3 Discrete analogue of HRSD distributions 

Steutel and van Harn (1979) had developed discrete SD (DSD) distributions. 
We now introduce RSD and HRSD for Z + -valued distributions. Some aspects 
of discrete HID (DHID) laws and generalized AR(1) models on Z + have been 
discussed in Satheesh et al. (20106). 

Definition 3.1 A PGF P(s) is DHID if for each p G (0, 1) there exists a PGF 
P p (s) such that 

P(S)= { a -(a-l)P p k (s)y/^ P= a' ^ 
Theorem 3.1 (Satheesh et al. (2010a)) A PGF P(s) is DHID iff 

P{S) = (1 - logW" (3 - 2) 
where k > integer and R(s) is a PGF that is DID. 

Definition 3.2 A PGF P(s) is discrete HRSD (DHRSD) if for each c G (0, 1] 
and p G [0, 1) there exists a PGF P c , p {s) such that 

P(s) = PcA*){P + (1 - P)P k (^ - c + cs)Y' k . (3.3) 
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Denoting the classes of DHRSD, DSD and DHID distributions by Cdhrsd, 
Cdsd and Cdhid we can proceed as in Section 2 to arrive at 

Proposition 3.1 We have Cdhrsd = Cdsd H Cdhid- Further, whenever 
P(s) G Cdhrsd-, the PGF P c , p (s) in (3.3) can be written as 

PcM = Pc(s).P P (l-c + cs) (3.4) 

where P p (s) and P c (s) are given by 

Pp(s) = {p + (i-l)P k (s)y/ k (3,6) 

Again, considering the Z + -valued analogue of the generalized AR(1) scheme 
(2.8) with 0, the binomial thinning operator in Steutel and van Harn (1979) 
we have the INAR(l) model 

Y = J €n ' h with P robabilit y P» ^ 3 ^ 

i=i ( Yli=i c + Z)f=i e n,ii with probability(l - p). 

Assuming stationarity we have the following Propositions as in Section 2. 

Proposition 3.2a If {l 7 ^} describes the stationary INAR(l) model (3.7) for 
each c G (0, 1] and p G [0, 1) then the distribution of {Y 1hi \ is DHRSD. 

Conversely, 

Proposition 3.3b If P Y (s) is a PGF that is DHRSD with P c , p (s) = P e (s) for 
each c G (0, 1] and p G [0, 1) then there exists a stationary INAR(l) model 
described by (3.7) with P Y (s) the PGF of {Y n ^} and P e {s) that of {e n>i }. 

Example 3.1 Negative binomial(^, A) distributions with PGF 

P{s) = (1 + A(1 1 _ a))1/fc , k > integer, A > (3.8) 

is DHID. Further since it is also DSD, this distribution is HRSD. 

Example 3.2 We may also proceed in a general frame work as done in Example 
2.2 to construct PGFs that are HRSD. 

Satheesh and Sandhya (2010) has proposed a further generalization of 
HRSD distributions based on the notion of MID distributions of Gnedenko 
and Korolev (1996). 
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